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Based on first-principles calculations, we systematically study the electronic, dielectric, and 
piasmonic properties of two-dimensional (2D) electride materials X 2 N (X=Ca, Sr). We show 
that both Ca 2 N and Sr 2 N are stable down to monolayer thickness. For thicknesses larger 
than 1-monolayer (1-ML), there are 2D anionic electron layers confined in the regions be¬ 
tween the [X 2 N]''' layers. These electron layers are strongly trapped and have weak coupling 
between each other. As a result, for the thickness dependence of many properties such as 
the surface energy, work function, and dielectric function, the most dramatic change oc¬ 
curs when going from 1-ML to 2-ML. For both bulk and few-layer Ca 2 N and Sr 2 N, the 
in-plane and out-of-plane real components of their dielectric functions have different signs in 
an extended frequency range covering the near infrared, indicating their potential applica- 
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tions as indefinite media. We find that bulk Ca2N and Sr2N could support surface plasmon 
modes in the near infrared range. Moreover, tightly-bounded plasmon modes could exist 
in their few-layer structures. These modes have significantly shorter wavelengths ( few tens 
of nanometers) compared with that of conventional noble metal materials, suggesting their 
great potential for plasmonic devices with much smaller dimensions. 

Electrides are a special kind of ionic solids with cavity-trapped electrons serving as the 
anions^’^. These electrons are spatially separated from the cations in a regular crystalline ar¬ 
ray, and are not bound to any particular atom, molecule, or bond. The first crystalline electride, 
Cs'''-(18-crown-6)2-e“, was synthesized in 1983^, and several other electrides have been success¬ 
fully discovered or predicted later on"^"*"^. The properties of electrides are closely connected to the 
topology and geometry of the cavities which confine the anionic electrons^. The early examples of 
electride demonstrate confinement of zero-dimensional cavities or one-dimensional weakly linked 
channels^. Recently, a new type of electride with two-dimensional (2D) confinement of anionic 
electrons was discovered in dicalcium nitride (Ca 2 N) which has a layered structure From for¬ 
mal valence consideration, each unit cell of Ca 2 N should have one excess electron. Experimental 
measurements on properties such as the transport parameters and work functions, combined with 
first-principles calculations, indeed proved that there is a built-in anionic electron layer confined 
between the calcium layers, which agrees well with the chemical formula [Ca 2 N]+ e“^^. This 
discovery of 2D electride materials (here 2D refers to the anionic electron confinement topology) 
has generated great research interest. First-principles calculations have predicted several other 2D 
electrides^^ '^, including the other two alkaline earth sub-nitrides Sr 2 N and Ba 2 N^^. Recent the- 
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oretical work also indicated that the monolayer Ca 2 N might be mechanically exfoliated from the 
bulk, meanwhile, its electron confined layer could be maintained, and suitable encapsulation layers 
were proposed to protect it in ambient environment^^. 

Motivated by these recent experimental and theoretical progress on 2D electrides, and also 
by the surge of research activities on 2D materials in recent years initiated by the discovery of 
graphene, in this work, we conduct a systematically investigation of the electronic, dielectric, 
and plasmonic properties of monolayer and few-layer alkaline earth sub-nitrides Ca 2 N and Sr 2 N. 
We find that besides Ca 2 N, Sr 2 N is also dynamically stable down to monolayer thickness. Their 
phonon spectra exhibit characteristic features of 2D materials. For thickness greater than 1-ML, 
besides the two surface electron bands, there are additional bands crossing the Fermi energy which 
are due to the anionic electron layers confined in the 2D interlayer regions between the [X 2 N]+ 
layers. These interlayer bands are lower in energy than the surface bands and are nearly degen¬ 
erate, indicating that they are only weakly coupled. The thickness dependence of properties such 
as surface energy, work function, and dielectric function are analyzed. The change in property 
with thickness is most dramatic when going from 1-ML to 2-ML, which is associated with the 
appearance of the first interlayer electron band. We find that Sr 2 N has a lower surface energy and 
lower work function compared with Ca 2 N. Due to their intrinsic structural anisotropy, these layered 
materials have highly anisotropic dielectric functions. In particular, the in-plane and out-of-plane 
components of dielectric function can have different signs, which occurs already for the bulk form 
in the near infrared frequency range. Moreover the dissipation is low in the range for the bulk ma¬ 
terials. This shows that bulk Ca 2 N and Sr 2 N could be ideal low-loss indefinite media^°. We further 
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investigate the possibility of guided surface plasmon modes for these metallic materials and show 
that bulk Ca 2 N and Sr 2 N have good performance comparable to the noble metals but operating at 
lower frequencies, while their thin film structure can have strongly bounded plasmon modes which 
have advantage of much shorter wavelength compared to the conventional noble metals. Our find¬ 
ings thus identify a class of promising 2D electride materials and reveal their potentials for future 
electronics and plasmonics applications. 

Results 

Crystal structure, dynamical stability, and surface energy. The alkaline earth subnitrides X 2 N 
(X=Ca, Sr, and Ba) in their bulk form can be synthesized by a direct solid-state reaction of the 
elements or through chemical reduction of the corresponding sesquinitrides^*. Their structures 
consist of alternating ABC-stacked (X-N-X) hexagonal layers, which are in the B2m space group 
(anti-CdCl 2 -type) with a high c/a ratio, where a and c are the in-plane and out-of-plane unit cell 
dimensions'^’Each (X-N-X) layer is closely packed. For example, the [Ca 2 N] layer has a 
thickness of 2.5lA which is much smaller than that of the ordinary fee Ca (111) layers (3.18A). 
This is typically attributed to the ionic bounding in the layer. Meanwhile, the separation between 
two neighboring layers is relatively large, leading to 2D confined spaces for the anionic electrons 

In this work, we investigate the properties of monolayer and few-layer thin films of these 
subnitrides (here one monolayer (ML) refers to one (X-N-X) unit layer). Since large interlayer 
separations (> S.SA) usually signal a possible role of long-range van der Waals interaction, the 
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van der Waals corrections was included in the first-principles calculation^"^. The calculation details 
are described in the Method section. As an example, the crystal structure of 3 -ML X2N (X=Ca, 
Sr) is schematically shown in Fig. la. The optimized lattice parameters for isolated Ca2N and Sr2N 
thin films with thickness from 1 -ML to 5 -ML are listed in Table 1 . One observes that despite slight 
variations with respect to the number of layers, the calculated in-plane lattice constants, layer 
thicknesses, and interlayer distances are close to the experimental values of the corresponding 
bulk structures (for bulk Ca2N, a = 3.62A, layer thickness= 2.5lA, and layer separation= 3.86A; 
for bulk Sr2N, a = 3.85A, layer thickness= 2.7lA, and layer separation= 4.19A)^^. Generally, 
the interlayer spacing is much larger than the layer thickness by about 50 %. Sr2N has a larger 
interlayer spacing (around 4A) than Ca2N due to the larger atomic number of Sr. 

The dynamical stability of 1 -ML Ca2N has been studied in Ref. 19 . Here we find that few- 
layer Ca2N and Sr2N (including Sr2N monolayer) are also dynamically stable through analysis of 
their phonon spectra. In contrast, the similar structures of 1 -ML and 2 -ML Ba2N are not stable 
due to the presence of imaginary frequencies in their phonon spectra. Therefore, we only focus on 
Ca2N and Sr2N in the present work. As representative examples, the phonon dispersions of mono- 
layer Ca2N and Sr2N are shown in Fig. 2 . The absence of imaginary frequencies in the whole Bril- 
louin zone demonstrates that the dynamical stability of the corresponding structures. The in-plane 
and out-of-plane transverse acoustic modes are not degenerate due to the structural anisotropy 
which is typical for layered materials. More importantly, while the in-plane transverse acoustic 
modes have a linearly dependence on their wavevectors in the vicinity of F-point, the out-of-plane 
acoustic (ZA) phonons exhibit a quadratic dispersion around F-point. The parabolic dispersion 
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of ZA modes is a characteristic feature of layered materials, which is consistent with the macro¬ 
scopic elastic theory of thin plates^^. This feature has been frequently observed in other layered 
2D materials, such as graphene^^, layered transition metal dichalcogenides^^, and phosphorene^*. 

We then calculate the surface energies of 1-ML to 5-ML Ca 2 N and Sr 2 N and the results 
are plotted in Fig.3. Because of their layered structure and the large interlayer separations, the 
surface energies are relatively low. The calculated values are comparable to the surface energy of 
graphene 12meV/A^)^^. Previous calculations on 1-ML Ca 2 N have suggested its possibility to 
be mechanically exfoliated from the bulk. From our result in Fig.3, the surface energy decreases 
with the increasing number of layers. Hence few-layer Ca 2 N would be more easily exfoliated 
compared with monolayer Ca 2 N, as naturally expected. Sr 2 N has an even larger interlayer spacing. 
Surface energy for 1-ML Sr 2 N is about 29meV/A^, and it decreases to about 20meV/A^ for 5-ML. 
With the same number of layers, Sr 2 N’s surface energy is smaller than that of Ca 2 N by about 
6 ~ SmeV/A^. This implies that Sr 2 N monolayer and few-layers could be more easily obtained by 
mechanical exfoliation using Scotch tapes or AFM tips. 

Electronic structure and 2 D confined electron layers. Bulk Ca 2 N and Sr 2 N have conduction 
electrons confined in the 2D interlayer regions playing the role of anions^^’^^. One key question 
is whether these anionic electron layers are maintained in their thin film structures. We find that 
the answer is positive. In Fig.4, we show the electronic band structures of X 2 N (X=Ca, Sr) from 
1-ML to 3-ML. First, one notices that all these systems are metallic, with Fermi level lying in 
partially filled dispersive energy bands. In 1-ML Ca 2 N (and Sr 2 N), there are two bands crossing 
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Fermi level (see Fig.4a and 4d). By inspecting the charge density distribution, one can verify that 
these two bands are from the 2D confined electron layers residing on the two sides of the [Ca 2 N]+ 
layer, and were referred to as the 2D electron gas in free space states^^. The two bands are energy 
splitted due to the coupling between the two sides. When a second [Ca 2 N] layer is added, i.e. 
for a 2-ML structure, there appears an additional band crossing Fermi level, which has similar 
dispersion as the two surface bands but is lower in energy (see Fig.4b). This band is from the 
confined anionic electron layer between the two [Ca 2 N]+ layers, as we will show in the following. 
The energy for this 2D band is lower because the confinement in the interlayer region is stronger 
than that for the surface, leading to a deeper potential well. One also observes that the splitting 
between the two surface bands is reduced, as a result of the reduced coupling between the two 
with increasing structural thickness. When going to 3-ML, there are two interlayer gap regions 
and indeed there is one more band crossing Fermi level (see Fig.4c). One observes that the two 
bands from the interlayer confined 2D electrons almost coincide, indicating that these states are 
strongly confined and coupling between neighboring interlayer regions is very small. This clearly 
demonstrates the 2D character of the confinement topology. The above features persist when more 
[Ca 2 N] layers are added. Then more interlayer confined anionic electron bands will appear and 
are nearly degenerate in energy. Their dispersion and bandwidth are almost independent on the 
thickness. For Sr 2 N, the band structure is quite similar to that of Ca 2 N with the same number 
of layers (with a slight decrease of the bandwidth, c.f. Fig.4d-f), implying that the 2D electride 
character is also maintained in its thin film form. 

To visualize the real space distribution of the electronic states, we take 2-ML Ca 2 N as an 
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example, and plot their partial electron densities for states within 0.05eV around Fermi energy 
Ep in Fig.5a. It is clear that the states around Fermi level are mainly located in three confined 
2D regions: outside the two surfaces and in the interlayer space between the [Ca 2 N]''' layers. We 
further plot the total charge densities of each of the three conduction bands for 2-ML Ca 2 N. As 
shown in Fig.5b-d, the lower band is from the 2D states in between the [Ca 2 N]+ layers, while the 
two higher bands are from the states confined at the surfaces. Similar result holds for Sr 2 N. This 
analysis validates our claims in the previous discussion. 

Electron-localization function (ELF) is useful for the analysis of the degree of electron lo¬ 
calization and the bonding character^*’’In Fig.6a and 6c, we show the ELF maps of 2-ML Ca 2 N 
and 2-ML Sr 2 N. Three delocalized electron layers can be seen. One observes that there is no bond¬ 
ing localization attractor between the confined electron layers and the [Ca 2 N]+ ([Sr 2 N]+) layers, 
indicating the bonding between them is of ionic type. When one valence electron is removed, the 
surface confined layers will be vacated while the interlayer anionic electron layer remain largely 
intact, as shown in Fig.6b and 6d. This feature also implies that the out-of-plane work function of 
these layered materials are mainly determined by the 2D electron layers on the surface. 

Previous studies found that bulk Ca 2 N has a highly anisotropic work function ranging from 
2.6eV (for (100) surface) to 3.4eV (for (001) surface)In Fig.7, we plot our results of the out- 
of-plane work function as a function of the number of layers for both Ca 2 N and Sr 2 N. One notes 
that there is a sharp decrease of the work function from 1-ML to 2-ML for both materials, due 
to the appearance of the strongly confined interlayer states that decouples the two surfaces. For 
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Ca 2 N, with increasing layer number, the work function approaches a value around 3.4eV, which 
is consistent with previous studies. The work function of Sr 2 N is less than that of Ca 2 N by about 
0.2eV, showing that the electrons in Sr 2 N are more loosely bound. We average the partial electron 
density around Fermi level in the layer plane. The resulting ID density profiles for 2-ML with 
and without one electron removed are plotted in Fig. 8. One observes that there are three peaks 
corresponding to the three anionic electron layers. After one electron is taken away, the two side 
peaks disappear while the central peak remains. The two side peaks of Sr 2 N have larger width than 
those of Ca 2 N, implying a higher degree of delocalization of these surface states hence a lower 
work function. 

Dielectric function. The first-principles calculation based on DFT has proved to be a powerful 
tool for the study of dielectric function for metals, including ultrathin metallic films down to few- 
layer thickness^^"^^. In crystalline solids, the dielectric function e{uj) consists of two contributions: 
a Drude-like intraband contribution and an interband contribution. The imaginary part of the inter¬ 
band contribution involves the interband matrix elements of the momentum operators, and can be 
evaluated directly in DFT^"^. Its real part can then be calculated via the Kramers-Kronig relation. 
The intraband contribution is typically treated by the Drude model, in which the plasma frequency 
(tensor) Up can be evaluated by DFT^"^. The calculation details are presented in the Method section. 

We first consider the dielectric properties of bulk Ca 2 N and Sr 2 N. Simple metals such as Au, 
Ag, and A1 generally have a large isotropic bulk plasma frequency around lOeV. In contrast, Ca 2 N 
and Sr 2 N are of layered structure, hence their dielectric properties are expected to exhibit intrinsic 


9 


anisotropy. We find that bulk Ca 2 N has an in-plane (i.e. in the layer plane) plasma frequency 
of 3.14eV and an out-of-plane plasma frequency cop^zz of 0.95eV; bulk Sr 2 N has an in-plane 
plasma frequency of 2.94eV and an out-of-plane plasma frequency of 1.05eV. These frequencies 
are significantly lower than those for simple metals, partly because the conducting electron density 
is lower. For each material, the out-of-plane plasma frequency is much lower compared with the 
in-plane plasma frequency, showing that the neighboring anionic electron layers are only weakly 
coupled. In Fig.9, we plot the real and imaginary part of in-plane and out-of-plane dielectric func¬ 
tions for the two materials. They show the character of metallic behavior, i.e. Drude peaks at 
low energy due to intraband contribution and the real part of e{u!) crossing from negative value to 
positive value with increasing frequency. The strong anisotropy in the dielectric function is also 
quite obvious. Remarkably, because of this anisotropy, the real part of in-plane dielectric function, 
Re£:a;a;(a;), changes sign at a frequency which is different from that for 'Reezz{oj), leading to an 
extended frequency window in which the two components have different signs (marked as the blue 
shaded regions in Fig.9). This sign difference is the characteristic feature of so-called indefinite 
media^*^, which was proposed in the study of metamaterials and has important potential applica¬ 
tions such as near-field focusing and building hyperlenses that can transform evanescent fields into 
propagating modes"^'*’"^*. Previous realizations of indefinite media are mostly in artificially assem¬ 
bled structures which require complicated fabrication process and usually have high dissipation. 
Our results suggest that crystalline solids Ca 2 N and Sr 2 N in their bulk form would just be indefinite 
materials for a frequency range spanning the near infrared. Moreover, one notes that the imaginary 
parts of and in this frequency range are very small (~ 0.1), implying that they could be 
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ideal low-loss indefinite materials. 


Next we focus on the dielectric functions of Ca 2 N and Sr 2 N few-layer structures. With 
decreasing film thickness, one expects that the anisotropy effects would be even stronger. The 
intraband contribution to the out-of-plane dielectric response becomes negligible, and it has 
been shown that for simple metals which are isotropic in bulk can exhibit strongly anisotropic e{(jj) 
down to few-layer thickness^^. In Fig. 10a, we plot the results of in-plane plasma frequency ojp^xx 
as a function of the number of layers for both materials. One observes that the frequency slowly 
decreases with increasing number of layers and approaches its bulk value. The plasma frequency 
of Sr 2 N is lower than that of Ca 2 N by about 0.2eV for each thickness. The square of plasma 
frequency is roughly proportional to the product of the carrier density and the inverse of effective 
mass. In Fig. 10b, one can see that from 1-ML to 2-ML, there is in fact a sharp decrease in the 
carrier density, as characterized by the DOS per unit volume around Fermi level, primarily due to 
the appearance of the large interlayer space of 2-ML. However, from the band structure in Fig.4b, 
one notices that the additional 2D electron band (marked in red color) has a smaller effective mass 
around Fermi level than the two surface bands, hence decreasing the average effective mass. This 
compensates the decrease in carrier density, resulting in an overall small change in the plasma 
frequency between 1-ML and 2-ML. 

In Fig.l 1, we show the interband contribution to the imaginary part of the in-plane dielectric 
function, which is closely connected to band structure features. One noticeable feature 

is the sharp peak around 0.3eV for both materials at 1-ML thickness. This peak has a sharp drop 
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from 1-ML to 2-ML and shifts to a lower energy. Its appearance can be attributed to the (almost) 
parallel sections of the two surface bands for 1-ML, as indicated in the red shaded regions in Fig.4a 
and 4d. From 1-ML to 2-ML, the splitting between two surface bands decreases and the volume 
increases a lot due to the appearance of the large interlayer region, resulting in the observed change 
of the peak. For thicknesses larger than 1-ML, there is another peak around leV, which is due to 
the transition between the interlayer 2D bands and the surface bands (c.f. Fig.4). From the band 
structures, it is clear that the response below leV is mainly contributed by the states from the 2D 
anionic electron layers. 

The total in-plane and out-of-plane dielectric function with their real and imaginary parts 
for Ca 2 N few-layers are plotted in Fig. 12. For Sxx^ the low energy part is dominated by the intra¬ 
band Drude-like contribution. Re£a:a:(a;) is negative below about leV. Meanwhile, the intraband 
contribution is negligible for ezz{<^)- Hence Re£ 22 (c(j) is largely positive for the low energy part. 
The sign difference between Re £^2 and 'Reexx in the low energy range again signals a possible 
indefinite material. This mechanism for the anisotropic dielectric function has been discussed in 
the case of Au ultrathin films^^. One observes that the change of e{uj) against thickness is most 
dramatic from 1-ML to 2-ML, due to the introduction of the first interlayer anionic electron layer. 
The imaginary part of dielectric function for 1-ML is much larger than that for 2-ML, particularly 
in the near infrared frequency range, which indicates more energy dissipation for possible plasmon 
modes as we discuss later. In addition, the variation of Sxx versus thickness is overall less dramatic 
compared with that of Szz, which again reflects the fact that the conducting electrons are strongly 
confined within the 2D in-plane layered regions. Due to the same reason that the different conduct- 
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ing 2D layers are only weakly coupled, the results for both in-plane and out-of-plane components 
converges rapidly with increasing thickness (note that the low energy part of the out-of-plane com¬ 
ponent would eventually be dominated by a Drude-like intraband contribution when approaching 
the bulk limit). One notes that there is only small difference between the results of 4-ML and 
5-ML. We also calculated 8-ML which is the largest system size within our computing capability. 
The obtained results are very close to those for 5-ML. Similar features are also demonstrated in 
the results for Sr 2 N few-layer structures, as shown in Fig. 13. 

Surface plasmon modes. Surface plasmon modes are confined electromagnetic excitations prop¬ 
agating at an interface between a conductor and a dielectric"^^. It typically requires a sign change of 
Ree across the interface. Therefore the most commonly used materials for plasmonic applications 
are metals such as Au, Ag, and A1 which have a range of frequencies with negative Re£:(a;)"‘^. 
Since Ca 2 N and Sr 2 N are also metals, one may also wonder whether they could also support sur¬ 
face plasmon modes. This question is particularly interesting when considering the frequency 
range in which the material shows indefinite medium property. In such case, one expects that the 
positive ReSzz component may compete with the negative component and tend to destroy 

the bounded plasmon modes. 

Let’s first consider the interface between bulk Ca 2 N or Sr 2 N and a dielectric medium charac¬ 
terized by a frequency-independent dielectric constant ej, > 0. The interface is parallel to the layer 
plane. As we have shown, Ca 2 N and Sr 2 N in the bulk form already have strong anisotropy in their 
dielectric functions. Assume that the interface supports a transverse magnetic (TM) plasmon mode 
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travelling along the interface with a wave vector jS. Following standard derivation using Maxwell’s 
equations"^^, one obtains that 



( 1 ) 


Here = w/c is the wave vector in vacuum, c is the speed of light, and Exx and Szz are the two 
components of the dielectric function for our conducting material. If Sxx = the result in Eq.(l) 


reduces to the familiar result for isotropic metal"^^. The dispersion characteristics of the surface 


plasmon modes for both materials are shown in Fig. 14, with sa = 2.25 (appropriate for Si02). 


In Fig. 14a, the modes lying to the right of the light line (in the dielectric medium) are bounded 


to the surface. One observes that the results show characteristic surface plasmon peaks similar 
to simple metals and confined plasmon modes still exist within the frequency range where the 
material shows indefinite medium property. A major difference is that for simple metals the peak 
position, known as surface plasmon frequency Wgp, occurs at higher energies, e.g. cOgp ~ 3.4eV for 
Ag/Si02 interface'^"^; while Wsp is much lower here, around l.leV-1.2eV in the near infrared range. 
The lmf5 shown in Fig. 14b is connected to the energy damping during the mode propagation. 
Since high dissipation occurs around Wsp, for practical applications, modes with frequencies less 
than ujsp are used. Here for Ca 2 N, if we take a; = 1.02eV, the corresponding surface plasmon 
wavelength Agp — 656nm, and we have a long propagation length L = l/(2Im/3) ~ 5.15/im and 
the decay length in the dielectric is Id = l/Re^/— Sdkp ~ 178nm. Longer propagation length 
can be achieved at lower frequencies, e.g. at a; = 0.64eV, we can have L ~ 173/im but the mode 
confinement is reduced, with Id ~ 737nm. These values are comparable to those for the noble 
metals (in the visible or UV frequency range) which are the usual building blocks for plasmonics 
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devices. However, the operating frequency here is much lower. From above discussion, we see 
that despite the intrinsic anisotropy, bulk Ca 2 N and Sr 2 N can still support surface plasmons and 
could be suitable plasmonics materials in the near infrared frequency range. 

We then turn to the thin films of Ca 2 N and Sr 2 N sandwiched between dielectric materials. 
Besides the change in the thickness-dependent dielectric function, an important effect in thin films 
is that the surface plasmon modes at two interfaces could couple and form two modes with opposite 
parity: a symmetric mode (L—) and an antisymmetric mode Their dispersions have been 

derived before in the study of ultrathin metallic films^^, and are quoted in our Method section. 
In Fig. 15, we plot the surface plasmon dispersion characteristics for Ca 2 N with different film 
thicknesses. One observes that the antisymmetric modes are lying on the light line of the dielectric 
medium, indicating that they are squeezed out of the metal region forming unbounded modes 
propagating in the dielectric material, which is similar to the case of Au ultrathin films. Meanwhile, 
pronounced plasmonic peaks do show up for the symmetric modes, clearly indicating that these 
are bounded surface plasmon modes. Again, one observes that the variation with thickness is most 
dramatic between 1-ML and 2-ML. For 1-ML, the bounded modes occur around leV. The high 
peak in Re/5 vs. to shows the modes are strongly bounded to the metallic layer. For thicknesses 
of 2-ML to 5-ML, the dispersions are quite close. The corresponding surface plasmon frequencies 
shift to around 1.2eV. As for the imaginary part of /5, 1-ML structure has relatively large values, 
whereas Im/5 for 2-ML to 5-ML almost collapse on a single curve and deceases rapidly in the range 
below 0.8eV. For Sr 2 N few-layers, the results for surface plasmon dispersion show similar features. 
The dispersions for symmetric (L—) modes are shown in Fig. 16. The antisymmetric (L+) modes 
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are again unbounded hence are not shown. 


Compared with their bulk results, one observes that the surface plasmon frequency cjgp is 
more or less the same, but both Re/3 and Im/3 are greatly increased by two orders of magnitude. 
This means that the plasmon wavelength and the confinement scale are much smaller, which are 
desired. However, the propagation length is also decreased at the same time. Therefore, it is 
more meaningful to consider the dimensionless ratio Re/3/Im/3, which measures how many sur¬ 
face plasmon wavelengths can be covered before the wave loses most of its energy, as well as the 
wave localization (or wave shrinkage) quantified by Aair/Agp, where Aair = 27rc/a; is the wave¬ 
length in air. In Fig. 17, we show these two dimensionless characteristics for the two materials 
as functions of Aair- The results for 3-ML and 4-ML are similar to those of 5-ML hence are not 
shown. One observes that for 1-ML Re/3/Im/3 is lower for most frequencies, due to the relatively 
high dissipation associated with Im/3 (c.f. Fig.lSd). Generally, the wave localization reaches its 
peak near the surface plasmon resonance. However, there Re/3/Im/3 is small due to the enhanced 
dissipation. While larger Re/3/Im/3 can be achieved for longer wavelengths, the wave localization 
becomes poor in that range. This tradeoff is a typical feature for surface plasmons. For application 
purposes, a compromise can be reached somewhere in between when the two are comparable. 

Since 1-ML Ca 2 N and Sr 2 N have considerably larger dissipation, in the following we take 
2-ML Ca 2 N and Sr 2 N as examples. For Aair = 2/im (0.64eV) in the near infrared range, 2-ML 
Ca 2 N has a decay length 1^ = 5.4nm in the dielectric, showing that the plasmon mode is strongly 
confined to the metallic layer. The ratios Aair/Agp = 56 and Re/3/Im/3 = 61. These two values are 
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considerably larger than what noble metals could achieve at their optimal working frequencies. The 
corresponding values for Sr 2 N at Aair = 2/rm (tc = 0.64eV) are given by Id = 4.4nm, Aair/Agp = 68 
and Re/3/Im/3 = 17. In comparison, the surface plasmon modes for noble metals such as Au 
and Ag are unbounded at this wavelength. Noticeably, the surface plasmon wavelengths in these 
materials are very small, Agp = 34.9nm for Ca 2 N and Agp = 28.9nm for Sr 2 N (at co = 0.64eV), 
while the shortest Agp for Au and Ag (occur in the visible light range) would be at least larger than 
lOOnm. These suggest that Ca 2 N and Sr 2 N thin films have great potential for making plasmonic 
devices operating in the near infrared range with much smaller scales. 

Discussion 

The loosely bound anionic electrons in electrides are highly reactive in ambient conditions. They 
could be used as good electron donors and as excellent catalysts for chemical reactions' ' . However, 
for physical applications, one needs to stabilize its property, e.g. by effective encapsulation. For 1- 
ML Ca 2 N, a possible encapsulation scheme using 2D insulating layers of graphane was proposed'^. 
In the case of few-layer alkaline earth subnitrides, since as we discussed each interlayer anionic 
electron layer are strongly confined in 2D regions and coupling between layers is small, one expects 
that the the most reactive electrons are from the surface layers, while the electron layers inside 
should be less reactive. A detailed study of this point and possible encapsulation schemes for Sr 2 N 
will be deferred to a future work. 

Compared with Ca 2 N, Sr 2 N has a smaller electrostatic potential associated with the larger 
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atomic number of Sr. This was reflected in the lower work function, and may also lead to a high 
electron mobility for Sr 2 N few-layers'^. Previous experimental studies have shown that bulk Ca 2 N 
has high mobility of 520cm^/(V s)^^. One expects that the mobility for Sr 2 N may be even higher. 
Hence Sr 2 N few-layers as new 2D conducting materials could have good potential for electronics 
applications. 

As having been demonstrated in other 2D layered materials, strain engineering has proved to 
be a powerful tool to modify and control the material properties. The layered materials discussed 
here have a large interlayer spacing (> 3.5A) between the [X 2 N]+ layers. We also expect that 
applying strain could be a good method to tune the material properties such as the dielectric func¬ 
tion, the plasmonic dispersion, and the carrier mobility. A systematic study of the strain effects is 
currently underway. 

Methods 

First-principles calculations. First-principles calculations were carried out using the Vienna ab- 
initio simulation package (VASP)"^^ "^®, based on the density functional theory (DFT). The exchange- 
correlation functional was treated using Perdew-Burke-Ernzerhof generalized gradient approximation"^’. 
The projector augmented wave (PAW) pseudopotential method"^^ is employed to model interactions 
between electrons and ions. The cutoff for plane-wave expansion is set to be 600 eV. The vertical 
distance between thin films (the thickness of the vacuum gap) is at least 18 A, which is large enough 
to avoid artificial interactions between the film and its periodic images. Both the atomic positions 
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and lattice constant were fully relaxed using conjugate gradient method. The convergence cri¬ 
teria for energy and force were set to be 10“^eV and O.OleV/A, respectively. DFT-D2 method 
was applied to describe the long-range van der Waals interaction. The Brillouin zone integrations 
have been carried out on a T-centered fc-mesh. Monkhorst-Pack k-point meshes"^^ with sizes of 
15x15x1 and 31x31x1 were used for geometry optimization and static electronic structure cal¬ 
culation, respectively. In the later study of optical properties, the sizes of fc-mesh are significantly 
increased to 61x61x1 and 41 x41 x9 for thin films and bulk respectively, to achieve highly con¬ 
verged results. For the integration over the Brillouin zone in calculating dielectric functions, we 
used the first order Methfessel-Paxton method^^ with a value of O.leV. The phonon dispersions of 
the structures were calculated by using density functional perturbation theory as implemented in 
PHONOPY code^O’^i. 


Calculation of dielectric functions. The optical properties of solids are mainly due to the re¬ 
sponse of the electron system to a time-dependent electromagnetic perturbation. For metals, the 
optical complex dielectric function consists of interband and Drude-like intraband contributions: 

£(CC) (CC). (2) 

The imaginary part of the interband part can be calculated using the results from DFT calculations 


Im (^)] ^ ^ y ‘^fnk{^m,k+qea\^n,k'){^n,k\^m,k+qep'} 


nm]k 


X 




( 3 ) 
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where a and /3 refer to Cartesian coordinates, ea{/ 3 ) are unit vectors, V is the volume of the unit cell, 
\un,k) and En,k are the periodic part of the Bloch wave function and the corresponding eigenenergy 
for band n and wave vector k, and fnk is the Fermi-Dirac distribution function. The real part 
of interband contribution can be obtained through the Kramers-Kronig relation. The intraband 
contribution is usually modeled by the Drude model: 


2 

T r^intra/, 


Re [4T(‘e)] = 1 - 


00 , 


p,ap 


(4) 

(5) 


a ;2 -|- 

Here 7 is a life-time broadening obtained either from a higher-order calculation or from exper¬ 
iments. In our calculation, we used the experimental determined electron life-time (of 0.6ps)'^ 
for bulk Ca 2 N to estimate 7 (~ l.lmeV). The same value was also used for the calculations of 
Sr 2 N due to their similar electronic structures. We have checked the sensitivity of our results’ de¬ 
pendence on 7 by repeating the calculations of e{oj) with 7 value varying from 7 /IO to 87 . The 
obtained results are of little difference. This is because the value of 7 is already quite small (re¬ 
flecting the fact that these materials are good metals). The oop^ap is the plasma frequency tensor 
which can be calculated using 


^p,al3 


47re^ 

~Y~ 





dEn^k \ 
dk ) 



( 6 ) 


The dielectric functions and plasma frequencies were suitably renormalized to exclude the vacuum 
region from the unit cell in our calculations. 


Surface plasmon modes calculations. For thin film structures, the plasmon modes at two sur¬ 
faces would couple and form two modes with different parity"^'^. The equations governing their 
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dispersions have been derived before^^: 


L+: exx\j- edkl - - s^J^tan = 0; (V) 

L-: exx\l- £dkl + - e^l cot = 0- (8) 

The first equation above is for the antisymmetric (L+) mode and the second equation is for the 
symmetric (T—) mode. Here L in the equations is the film thickness. Other quantities in these two 
equations are defined in the main text. In our calculation, we solve the two equations numerically 
using a two-dimensional unconstrained Nelder-Mead minimization algorithm^^ with a tolerance of 
10 “^^nm“^ in the complex wave vectors. 
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Table 1: Lattice parameters of Ca 2 N and Sr 2 N thin films with thicknesses from 1-ML 
to 5-ML. Here a is the in-plane lattice constant, t is the thickness of the film (vertical 
distance between the top and bottom atomic layers), L(1,2,3) refer to the thicknesses of 
the [Ca 2 N]+ layers, L1 Is for the outermost layer(s), L2 is for the next outermost layer(s) 
and so on (see Fig.1), and Gij refers to the thickness of the interlayer gap region between 
Layer i and Layer j. 


No. of Layers 

a (A) 

t(A) 

LI (A) 

L2 (A) 

L3 (A) 

G12 (A) 

G23 (A) 

Ca2N 








1 

3.562 

2.516 

2.516 





2 

3.562 

8.608 

2.509 



3.598 


3 

3.552 

14.957 

2.514 

2.504 


3.712 


4 

3.550 

20.099 

2.516 

2.507 


3.721 

3.793 

5 

3.555 

26.419 

2.513 

2.505 

2.508 

3.703 

3.823 

Sr2N 








1 

3.778 

2.761 

2.761 





2 

3.775 

9.427 

2.757 



3.913 


3 

3.770 

16.260 

2.760 

2.748 


3.983 


4 

3.765 

21.821 

2.764 

2.754 


4.016 

4.068 

5 

3.768 

28.562 

2.758 

2.751 

2.753 

3.989 

4.061 
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Figure 1 a, Schematic figure of crystal structure of 3-ML X 2 N (X=Ca, Sr) viewed from 
[110] direction. It consists of three (X-N-X) unit layers with ABC-stacking. The symbols 
LI, L2, and G12 refer to the thicknesses of outermost layer, the next outermost layer, and 
the interlayer gap between LI and L2, respectively, as explained in the caption of label 1 . 
b, 2D Brillouin zone of X 2 N (X=Ca, Sr) few-layer structures in which the high symmetry 
points are labeled. 


Figure 2 Phonon dispersions, a. Phonon dispersion of 1-ML Ca 2 N. b. Phonon disper¬ 
sion of 1-ML Sr 2 N. 


Figure 3 Surface energies for few-layer Ca 2 N and Sr 2 N as a function of thickness d from 
1-ML to 5-ML. 


Figure 4 Eiectronic band structures of few-iayer Ca 2 N and Sr 2 N. a-c. Band structure 
of Ca 2 N with a 1-ML, b 2-ML, and c, 3-ML thickness, d-f, Band structure of Sr 2 N with d 
1-ML, e 2-ML, and f, 3-ML thickness. Fermi energy is set at zero. The two green colored 
bands are mainly from the 2D electron layers confined to the surface. The red colored 
bands (for films thicker than 1-ML) are mainly from the 2D electron layers confined in the 
interlayer regions. The red shaded rectangles in a and d indicate the regions contribute to 
the peaks ~0.3eV in interband Im£(a;) for 1-ML Ca 2 N and Sr 2 N, as observed in Fig4.(a) 
and Fig.(d). 
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Figure 5 Electron density distribution in 2-ML Ca 2 N. a, Partial electron density iso¬ 
surfaces (with value of 0.0003/Bohr^) for states in the energy range \E - Ef\ <0.05eV 
shown for a conventional unit cell, b-d, Band-decomposed electron density isosurfaces 
(with value of 0.003/Bohr^) for the three (colored) bands which cross Fermi level as shown 
in Fig.4(b). b is for the upper green band, c is for the lower green band, and d is for the 
red band. 

Figure 6 Electron localization function (ELF) maps for 2-ML Ca 2 N and Sr 2 N. a and 

c are the ELF maps for Ca 2 N and Sr 2 N respectively, shown for a conventional unit cell, b 
and d are for [Ca 2 N]+ and [Sr 2 N]+ respectively, where one valence electron is removed. 

Figure 7 Work functions for few-layer Ca 2 N and Sr 2 N as a function of thickness d from 
1-MLto 5-ML. 

Figure 8 Partial density of states for the energy range \E - Ef\ <0.05eV averaged 
over the ab plane for a 2-ML Ca 2 N and b 2-ML Sr 2 N. In each figure, the blue curve and 
the red curve are for the distributions before and after one valence electron is removed, 
respectively. The green shaded regions Indicate the locations of the [Ca 2 N]+ or [Sr 2 N]+ 
layers. 

Figure 9 Anisotropic dielectric functions for bulk X 2 N: a for Ca 2 N and b for Sr 2 N. 
The real and imaginary parts of in-plane component and out-of-plane component 
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are plotted using different colors. The blue shaded region in each figure indicates the 
frequency range in which Ree^^ and Ree^^ have different signs. 

Figure 10 a, In-plane plasma frequency for Ca 2 N and Sr 2 N thin films as a function 
of thickness d from 1-ML to 5-ML. b, eiectron density of states (DOS) of Ca 2 N and Sr 2 N 
at Fermi energy as a function of film thickness. In each figure, the dashed lines indicate 
the corresponding vaiues for the bulk. 

Figure 11 The interband contribution to for a Ca 2 N and b Sr 2 N thin fiim struc¬ 

tures from 1-ML to 5-ML. Resuits for different thicknesses are piotted using different col¬ 
ors. 

Figure 1 2 Dieiectric functions for Ca 2 N few-iayers. a and b are for the real and imag¬ 
inary parts of in-piane component respectively, c and d are for the real and imaginary 
parts of out-of-piane component respectiveiy. Results for different thicknesses are 
piotted using different colors. 

Figure 13 Dielectric functions for Sr 2 N few-iayers. a and b are for the real and imag¬ 
inary parts of in-piane component e^x respectively, c and d are for the real and imaginary 
parts of out-of-piane component respectiveiy. Results for different thicknesses are 
piotted using different colors. 
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Figure 14 Dispersion reiation of surface plasmon modes for an interface between a 
dielectric medium with = 2.25 and bulk X 2 N (X=Ca, Sr). The light line in the dielectric 
medium is also shown in grey coior. 


Figure 15 Dispersion characteristics of surface piasmon modes for Ca 2 N few-iayers in 
a dieiectric medium with sa = 2.25. a and b are for the antisymmetric mode (L+ mode), 
c and d are for the symmetric mode (L- mode), a and c show the reai part of the wave 
number component, b and d show the imaginary part of the wave number component. 


Figure 16 Dispersion characteristics of L- symmetric surface plasmon modes for Sr 2 N 
few-layers in a dieiectric medium with ea = 2.25. a shows the real part of the wave number 
component, b shows the imaginary part of the wave number component. 


Figure 17 Characteristics RejS/lmlS and Aair/Asp for the symmetric surface plasmon 
modes in few-layer a Ca 2 N and b Sr 2 N, plotted as functions of the waveiength in air Aair- 
Asp is the surface plasmon wavelength. The dielectric medium is with = 2.25. 
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